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^ : Abstract 

' We address several estimation problems in quantum optics by means of the 

maximum-likelihood principle. We consider Gaussian state estimation and the 
determination of the coupling parameters of quadratic Hamiltonians. More- 
over, we analyze different schemes of phase-shift estimation. Finally, the 
^ ' absolute estimation of the quantum efficiency of both linear and avalanche 

^ ' photodetectors is studied. In all the considered applications, the Gaussian 

. bound on statistical errors is attained with a few thousand data. 

o 

g : I. INTRODUCTION 



p • In order to gain information about a physical quantity one should, in principle, measure 

the corresponding quantum observable. In cases when the measurement can be directly 
implemented the statistics of the outcomes is governed (in ideal conditions, i.e. neglecting 
^ \ thermal, mechanical or other sources of classical noise) only by the intrinsic fluctuations of 
^1 the observable, namely by the quantum nature of the system under investigation. In practice, 
however, it is most likely that the desired observable does not correspond to a feasible 
^ . measurement scheme, or the physical quantity does not correspond to any observable at all. 

In such case one has to infer the value of the quantity of interest from the measurement 
of a different observable, or generally of a set of observables. In this situation, even in 
ideal conditions, the indirect parameter estimation gives an additional uncertainty for the 
estimated value, and the quantum estimation theory provides a general framework to 
optimize the inference procedure. 

In the recent years, the indirect reconstruction of observables and quantum states has 
received much attention. Among the many reconstruction techniques, the most successful 
is quantum homodyne tomography [0], which, indeed, is the only method which has been 
experimentally implemented [Q. Quantum tomography provides the complete characteri- 
zation of the state, i.e. the reconstruction of any quantity of interest by simple averages 
over experimental data. In many cases, however, one may be interested not in the complete 
characterization of the state, but only in some specific feature, like the phase or the ampli- 
tude of the field. Moreover, one can address the problem of characterizing an optical device, 
rather than a quantum state, like measuring the coupling constant of an active medium or 
the quantum efficiency of a photodetector. In all these cases, the desired parameter does 
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not correspond to a measurable observable, and contains only a partial information about 
the quantum state of light involved in the process. Our goal is to link the estimation of such 
parameters with the results from feasible measurement schemes, as homodyne, heterodyne 
or direct detection, and to make the estimation procedure the most efficient. 

Among all possible procedures for parameter estimation, the maximum-likelihood (ML) 
method is, in the sense discussed below, the most general, and widely usable in practice. The 
ML procedure answers to the following question: which values of the parameters are most 
likely to produce the results which we actually observe in the measurement ? This statement 
can be quantified, and the resulting procedure is the ML estimation of the parameters. 

Recently, the ML principle has been applied to the reconstruction of the whole state 
of a generic quantum system 1^,^. In that case the parameters of interest are the matrix 
elements of the density operator in a suitable representation. Bayesian and ML approaches 
have been also applied in neutron interferometry [^. 

In this paper, we focus our attention on the determination of specific parameters which 
are relevant in quantum optics, and analyze their ML estimation procedure in some details. 

In the next Section we briefiy review the ML estimation procedure as well as the method 
to evaluate its precision. In Section |II| we consider the estimation of the parameters of a 
Gaussian state and of the coupling constants of a generic quadratic single-mode Hamiltonian. 
As we will show, the two estimation problems are closely related, and ML principle leads 
to a fully general solution. In Section |V| we study different schemes of phase estimation, 
whereas in Section the ML principle is applied to the estimation of the quantum efficiency 
of both linear and avalanche photodetectors. Section ^ closes the paper by summarizing 
our results. 



II. MAXIMUM-LIKELIHOOD ESTIMATION 

Here we briefiy review the theory of the maximum-likelihood (ML) estimation of a single 
parameter. The generalization to several parameters is straightforward. Let p{x\X) the 
probability density of a random variable x, conditioned to the value of the parameter A. 
The form of p is known, but the true value of the parameter A is unknown, and will be 
estimated from the result of a measurement of x. Let Xi,X2, ■■■,xp^ be a random sample of 
size N. The joint probability density of the independent random variable Xi,X2, ■■■,xn (the 
global probability of the sample) is given by 

C{xi,X2, ...,XAr|A) = Il^^-^p{Xk\X) , (1) 

and is called the likelihood function of the given data sample (hereafter we will suppress the 
dependence of £ on the data). The maximum-likelihood estimator Aml = Kihiixk}) (MLE) 
of the parameter A maximizes C{X) for variations of A, namely it is given by the solution of 
the equations 

^ = 0; ^<0. (2) 
Since the likelihood function is positive the first equation is equivalent to 
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where 



N 



L{X)=\ogC{X) = Y.\ogp{xk\X) 



(4) 



k=l 



is the so-called log-likelihood function. The form of the ML principle in Eq. is the most 
often used in practice. 

The importance of MLE stems from the following theorems 

1. Maximum-likelihood estimators are consistent, i.e. they converge in probability to the 
true value of the parameter for increasing size of the data sample. 

2. The distribution of MLEs tends to the normal distribution in the limit of large samples, 
and MLEs have minimum variance. For finite samples the variance is governed by the 
Cramer- Rao bound (see below). 

There are also situations in which the MLE gives a poor estimation for a parameter. How- 
ever, for the distributions considered here the ML procedure is statistically efficient. 

In order to obtain a measure for the confidence interval in the determination of Aml we 
consider the variance 
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Y[dxkp{xk\\) 



L k 



(5) 



Upon defining the Fisher information 



it is easy to prove |Ty] that 



dx 



dp{x\\) 
dX 



p{x\X) 



(6) 



N F 



(7) 



where N is the number of measurements. The inequality in Eq. is known as the 
Cramer-Rao bound P| on the precision of ML estimation. Notice that this bound holds 
for any functional form of the probability distribution p{x\X), provided that the Fisher 
information exists VA and d\p{x\X) exists Vx. When an experiment has "good statistics" 
(i.e. a data sample large enough) the Cramer-Rao bound is saturated. As we will show in 
the following, the application of the ML principle in quantum optics generally corresponds 
to estimators for which the Cramer-Rao bound is attained with a relatively small number 
of measurements, i.e. the ML procedure provides an efficient estimation of the parameters. 
In Sections |ITl| and |I^ examples will be examined where the probability p{x\X) is Gaussian 



versus x and not Gaussian versus the parameter A, whereas in Section an example with 
discrete measurement outcomes {x = 0, 1) will be also analyzed. 
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III. GAUSSIAN-STATE ESTIMATION 



In this section we apply the ML method to estimate the quantum state of a single-mode 
radiation field that is characterized by a Gaussian Wigner function. Such kind of states 
represents the wide class of coherent, squeezed and thermal states. Apart from an irrelevant 
phase, we consider the Wigner function of the form 



W{x,y) 



2A2 



TT 



exp {-2A2 [e-^'-{x - Refif + e^'{y - Im/i)^] } , 



(8) 



and we apply the ML technique with homodyne detection to estimate the four real parame- 
ters A, r, Re/i and Im/x. The four parameters provide the number of thermal, squeezing and 
coherent-signal photons in the quantum state as follows 







nth = 


-(-- 

2 VA2 




sinh^ r , 


ncoh ~ 





(9) 

In terms of density matrix, the state corresponding to the Wigner function in Eq. (H) writes 



g = D{^) S{r) 



1 



nth + 1 \nth + 1 



nth 



S\r)D\ii) 



(10) 



where S{r) = exp[r(a^ — a^^)/2] and -D(/i) = exp{fia^ — fi*a) denote the squeezing and 
displacement operators, respectively. 

The theoretical homodyne probability distribution at phase with respect to the local 
oscillator is given by the Gaussian |[TT[] 



2A2 



vrfe^'' cos^ -|- e sin^ - 



X exp 



2A2 



e^** cos^ 6 + e 2'" sin^ - 



X — Re(/i e 



(11) 



From Eqs. (^) and (IT^) one easily evaluates the log-likelihood function for a set of 
homodyne outcomes Xi at random phase (pi as follows 



N 



2A2 



i=l 



TC[e^' COS 
2A2 



e^^' cos^ : 



e sin^ - 



Xi-Re(/ie-^'^') 



(12) 



The ML estimators A^l, ''"ml, Re/iML and Im/i^L are found upon maximizing Eq. (|T^) versus 
A, r, Re/i and Im/i. 

In order to obtain a global estimation of the goodness of the state reconstruction, we 
evaluated the normalized overlap O between the theoretical and the estimated state 
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0= ■ (13) 

/Tr[f.2]Tr[p2 i 



-MLJ 



Notice that O = 1 q = ^^ml- Through some Monte-Carlo simulations, we always found 
a value around unity, typically with statistical fluctuations over the third digit, for number 
of data samples = 50000, quantum efficiency at homodyne detectors rj = 80%, and state 
parameters with the following ranges: rifh < 3, ricoh < 5, and n^g < 3. Also with such 
a small number of data samples, the quality of the state reconstruction is so good that 
other physical quantities that are theoretically evaluated from the experimental values of 
Aml, ''"ml, Re/iML and Im/^ML are inferred very precisely. For example, we evaluated the photon 
number probability of a squeezed thermal state, which is given by the integral 

with C(0, nth, r) = {nth + sin^ (j) + e^*" cos^ 0) + |. The comparison of the theoretical 

and the experimental results for a state with nth = 0.1 and ngq = 3 is reported in Fig. |l|. 
The statistical error of the reconstructed number probability affects the third decimal digit, 
and is not visible on the scale of the plot. 

The estimation of parameters of Gaussian Wigner functions through the ML method 
allows one to estimate the parameters in quadratic Hamiltonians of the generic form 

H = aa + a*a^ + ipa^a + ^^a^ + ^C(^^^ ■ (15) 

In fact, the unitary evolution operator U = e~*^* preserves the Gaussian form of an input 
state with Gaussian Wigner function. In other words, one can use a Gaussian state to probe 
and characterize an optical device described by a Hamiltonian as in Eq. (p!5|). Assuming 
t = 1 without loss of generality, the Heisenberg evolution of the radiation mode a is given 
by 

all = -fa + Sa^ + jj. , (16) 

with 



7 = cos(v'^2 _ |^|2) _ ^—^= sin(v'^2 _ |^|2) 
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For an input state g with known Wigner function Wg{P,j3*), the corresponding output 
Wigner function writes 

Wu,u^ {(3 , (3*) = W,[{(3 - /i)7* - {[3* - /i*)5 , {f3* - /i*)7 - {f3 - /i)5*] . (18) 

Hence, by estimating the parameters 'y,S,fi and inverting Eqs. (0), one obtains the ML 
values for a,(p, and ^ of the Hamiltonian in Eq. (|TB|). The present example can be used 



in practical applications for the estimation of the gain of a phase-sensitive amplifier or 
equivalently to estimate a squeezing parameter. 
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IV. PHASE ESTIMATION 



The quantum-mechanical measurement of the phase of the radiation field is the essential 
problem of high sensitive interferometry, and has received much attention in quantum optics 
jr^ . The problem arises because for a single mode of the electromagnetic field there is no 
selfadjoint operator for the phase, hence a more general description of the phase measurement 
is needed on the ground of estimation theory 

In the following we apply the ML method to different schemes of phase estimation and 
evaluate the corresponding sensitivity. 

A. Heterodyne detection on coherent state 

For a coherent state with amplitude Ae^'^ the probability density for complex outcome 
at the i-th heterodyne measurement is given by 

= -exp(-|ai - Ae^^n . (19) 
vr 

The max-likelihood condition dL/dip = provides the MLE for the phase ip. One obtains 
V'ml = arg(a), where the overline denotes the experimental average over N heterodyne 
outcomes, namely a = (J2f=i o^'i) / N . For small phase-shift -i/^ ~ the Cramer-Rao bound 
gives the constraint > l/y/2nN, n being the average photon number (n = A^). 

B. Homodyne detection at random phase on coherent state 

In this case the homodyne probability for outcome Xi at the i-th measurement at phase (pi 
writes 



p(xi,0i) = ^exp{-2[xi-Acos(0i-V^)]2} . (20) 

The ML condition provides for the estimator of ip the solution 

-f/^ML = arctg (^a;sin(/)/a;cos0) . (21) 

Also in this kind of phase-detection strategy, the variance of the estimator for small phase- 
shift satisfies 
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C. Homodyne detection at fixed phase on squeezed states 



The use of squeezed states and homodyne detection at the phase corresponding to the 
squeezed quadrature offer a better resuh in terms of sensitivity. Consider the problem 
of estimating the phase ip in the state D{Ae'^'^)S{r)\0) with A,r > 0. The homodyne 
probabihty of outcome yi for the measurement of the quadrature Y = {a — a^)/2i writes 



I2e- 



2r 



TT 



exp 



(23) 



The MLE for ip is then given by ipML = arcsin(?//A). For small phase shift ~ the 
Cramer-Rao bound provides the relation 



2r ' 



(24) 



Upon maximizing the product A e versus the total number of photons in the state n 
+ sinh r, one obtains the optimal squeezing 



.2r 




(25) 



Notice that for A > 1, Eq. 



requires that an equal number of squeezing and coherent 



photons contributes to the total average power in the radiation, namely rtcoh — — 



In this case Eq. (|2^ ) rewrites 



9 1 



n/2. 



(26) 



namely one obtains the ideal limit for the sensitivity of phase estimation [jl],^. The bounds 
on sensitivity obtained in the previous examples are saturated within a rather small number 
of data samples. In Fig. |] we compare the experimental error obtained by a Monte-Carlo 
simulation of homodyne detection on squeezed states using 5000 data samples with the 
theoretical bound of Eq. (|^). We fixed the total number of photons at the value n = 50, 
and varied the squeezing fraction nsq/n. Notice how experimental and theoretical data 
compare very well. We estimated the statistical errors in Figs. 2-4 from the raw data by 
propagation of the errors on the evaluation of y, namely 



dy 



(27) 



Notice that for large data samples, ay —>■ e^'^''' /AN, and one recovers Eq. (^4]). As shown in 
Figs. 2-4, our estimation of errors approaches the Cramer-Rao bound, hence proving that 
the ML method for the phase estimation is statistically efficient. At the optimal value of 
squeezing fraction [see Eq. (|^)], the behavior oc 1/n is well reproduced, also at the small 
number 5000 of data samples, as shown in Fig. ^ 

Unfortunately, the result in Eq. (^) is very sensible to the effect of less-than-unity 
quantum efficiency 1] of realistic homodyne detectors. For 1] < 1, the homodyne probability 
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is given by a convolution of the ideal distribution in Eq. ( p3D with a Gaussian with variance 
(1 — ri)/4:ri. In such case, Eq. (|2^ ) is replaced by 



a. 



The optimal value of the squeezing factor e to minimize Eq. ( pS]) at fixed total number 
of photons is given by the solution in the interval / = [0, 1] of the cubic equation 

a;3 +(4A^ + x^-x- = . (29) 

V V J V 

Compare Fig. ^ with Fig. ^ where quantum efficiency t] = 0.8 has been used. Indeed, 
the optimal squeezing fraction rapidly approaches zero as r ~ rj for rj ^ 0. Such a 
detrimental effect of quantum efficiency is similar to the effect of losses in squeezed-state 
homodyne communication channels [|T3]. However, it can be partially stemmed by adopting 



a feedback-assisted homodyne detection 



V. ABSOLUTE ESTIMATION OF THE QUANTUM EFFICIENCY 

In principle, in a photodetector each photon ionizes a single atom, and the resulting 
charge is amplified to produce a measurable pulse. In practice, however, available photode- 
tectors are usually characterized by a quantum efficiency lower than unity, which means 
that only a fraction of the incoming photons lead to an electric pulse, and ultimately to 
a "count". If the resulting current is proportional to the incoming photon flux we have a 
linear photodetector. This is , for example, the case of the high flux photodetectors used in 
homodyne detection. On the other hand, photodetectors operating at very low intensities 
resort to avalanche process in order to transform a single ionization event into a recordable 
pulse. This implies that one cannot discriminate between a single photon or many photons 
as the outcomes from such detectors are either a "click", corresponding to any number of 
photons, or "nothing" which means that no photons have been revealed. In this section we 
apply the ML principle to the absolute estimation of the quantum efficiency of both linear 
and avalanche photodetectors. We suppose to have at our disposal a known reference state 
and, from the results of a measurement upon such a state, we infer the value of the quantum 
efficiency. 

Let us first study the case of linear photodetectors. As a reference state we consider a 
squeezed-coherent state, measured by homodyne detection. The effect of nonunit quantum 
efficiency rj on the probability distribution of homodyne detection is twofold. We have both 
a rescaling of the mean value and a broadening of the distribution. For a squeezed state 
|xo, r) = D{xo)S{r)\0) with the direction of squeezing parallel to the signal phase and to the 
phase of the homodyne detection (without loss of generality we set this phase equal to zero 



and xo, r > ) we have 15 



pJx) = , exp 



(x — 77x0)" 



A' 



2 



2A2 

\ (e-^^ + . (30) 
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The total number of photons of the state is given by n = Xq + sinh^ r, whereas the squeezing 
fraction is defined as 7 = sinh^r/ra. Apart from an irrelevant constant, the log-likelihood 
function can be written as 



L{r]) = log A'^ + ^ (x'^ + -qxl - 2r]XoX^ . (31) 



The resulting MLE is thus given by 




r^ML = 1 + e"'" + ^\l-Jl + 6Axl \x^ + (1 + e-^^){xo - 2x + xoe-2'^)a;ol • (32) 



A set of Monte Carlo simulated experiments confirmed that the Cramer-Rao bound is at- 
tained. The performances of the ML estimation can be compared to the "naive" estimation 
based only on the measurement of the mean value, i.e. t^av = x/xq. We expect this method 
to be less efficient, since the quantum efficiency not only rescales the mean value, but also 
spreads the variance of the homodyne distribution in Eq. (|5DD. In Fig. ^, on the basis of a 
Monte Carlo simulated experiment, we compare the ML and the average- value methods in 
estimating the quantum efficiency through homodyne detection on a squeezed state. The 
advantages of ML method are apparent, especially for the estimation of low values of t]. 
On the other hand, for small values of the squeezing fraction the two methods have similar 
performances, except for very low signals, whereas the ML estimation performs better. 

Let us now consider avalanche photodetectors, which perform the ON/OFF measurement 
described by the two-value POM 

00 

HoFF = E(l - \P){P\ Hon = I - Hopp , (33) 

p=0 

where I denotes the identity operator. Indeed, for high quantum efficiency (close to unity) 
IIoFF and Hon approach the projection operator onto the vacuum state and its orthogonal 
subspace, respectively. With avalanche photodetectors we have only two possible outcomes: 
"click" or "no clicks" which we denote by "1" and "0" respectively. The log-likelihood 
function is given by 

L{r]) = {N- N,) logPo(r?) + A^clog[l - Po{v)] , (34) 

where Po{r]) = Tr[^?noFF] is the probability of having no clicks for the reference state described 
by the density matrix g, N is the total number of measurements, and Nc is the number of 
events leading to a click. The maximum of L{ri), i.e. the MLE for the quantum efficiency, 
satisfies the equation 

PoiVu.) = 1 - ^ ' (35) 

whose solution, of course, depends on the choice of the reference state. The optimal choice 
would be using single-photon states as a reference. In this case, we have the trivial result 
'7ml = Nc/N. However, single-photon state are not easy to prepare [|l6l, and generally one 
would like to test r] for coherent pulses \a). In this case, we have Po(^) = exp(— japr]) and 
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?7ml = -rn^log fl - ^ ) . (36) 



The Fisher information is given by 



and therefore, for a weak coherent reference one has 



' ^^^"^^ , (37) 



er?|ap _ 1 

and 



- ^ (38) 



VI. SUMMARY AND CONCLUSIONS 

In quantum optics, there are several parameters of great interest corresponding to quanti- 
ties that are not directly observable. Among these, we studied the parameters of a Gaussian 
state, the phase of a squeezed-coherent state, and the quantum efficiency of either linear 
or single-photon resolving photodetector. In this paper, we have applied the maximum- 
likelihood method to the determination of these parameters using feasible detection schemes. 
In particular, we have considered homodyne detection and ON / OFF photodetection. In all 
cases here analyzed, the resulting estimators are efficient, unbiased and consistent, thus 
providing a statistically reliable determination of the parameters of interest. Moreover, by 
using the ML method only few thousand data are required for the precise determination 
of parameters. We stress that the ML procedure used in this paper can be applied to a 
broad class of estimation process, since it applies to any probability distribution p{x\X), as 
long as its functional form is known and the maximum of the likelihood function is unique. 
In conclusion, for the measurement of parameters pertaining to quantum states or optical 
devices, the ML procedure should be taken into account, in order to optimize data analysis 
and thus reducing the experimental efforts. 
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FIG. 1. Photon- number probability of a 
squeezed-thermal state. The black histogram for 
the theoretical, the gray one reconstructed state 
by means of the maximum likelihood method 
and homodyne detection. Number of data sam- 
ples N = 50000, quantum efficiency rj = 80%, 
number of thermal photons n^/j = 0.1, number 
of squeezing photons Ugq = 3. The statistical er- 
ror affects the third decimal digit, and it is not 
visible in the scale of the plot. 
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FIG. 2. Estimation of the statistical er- 
ror (circles) for the phase-shift measurement 
through the maximum likelihood method on a 
squeezed state of radiation, for different values 
of the degree of squeezing. The total number 
of photon of the state is fixed at n = 50. The 
solid line represents the Cramer-Rao bound on 
the errors. Only 5000 homodyne data have been 
used, and the bound is saturated, thus proving 
the efficiency of the method. 
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FIG. 3. Phase sensitivity vs total number of 
photons achievable through homodyne detection 
on squeezed states and maximum likelihood esti- 
mation, with optimal fraction of squeezing pho- 
tons [see Eq. (p5|)]. Compare the results of 
a Monte-Carlo simulation with 5000 homodyne 
outcomes (circles) with the theoretical behavior 
(solid line). 
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FIG. 5. Estimation of the quantum efficiency 
of linear photodetectors through homodyne de- 
tection on a squeezed state. Both plots report 
the ratio between the estimated value of the 
quantum efficiency and the true value, as a func- 
tion of the true value, (a) results obtained us- 
ing the maximum-likelihood method; (b) results 
by the "naive" average- value method. The ho- 
modyne sample consists of 50 blocks of 50 data 
each, whereas the reference state is a squeezed 
state with mean number of photons n = 1 and 
squeezing fraction 7 = 99% (nearly a squeezed 
vacuum) . 
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